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ABSTRACT 

In order to reconstruct the initial conditions of the universe it is important to devise a method that 
can efficiently constrain the shape of the power spectrum of primordial matter density fluctuations in 
a model-independent way from data. In an earlier paper we proposed a method based on the wavelet 
expansion of the primordial power spectrum. The advantage of this method is that the orthogonality 
and multiresolution properties of wavelet basis functions enable information regarding the shape of 
Pin(fc) to be encoded in a small number of non-zero coefficients. Any deviation from scale-invariance 
can then be easily picked out. Here we apply this method to simulated data to demonstrate that it 
can accurately reconstruct an input Pin(fc), and present a prescription for how this method should be 
used on future data. 

Subject headings: Cosmology; Cosmic Microwave Background, methods: data analysis 



1. INTRODUCTION 

Cosmic microwave background (CMB) and large scale 
structure (LSS) data are powerful probes of the universe. 
They are sensitive not only to the values of cosmological 
parameters but also to the initial conditions of the uni- 
verse that may have arisen due to inflation. While the 
general concept of inflation solves many problems in cos- 
mology and provides a framework within which the data 
are very well fit, not much is known about the details of 
inflation. 

The power spectrum of density perturbations that 
results from inflation is usually assumed to be scale- 
invariant, or sometimes a power-law or a power-law with 
a running tilt, and such parametrizations flt the data 
well so far. Simple models of inflation generically predict 
Gaussian adiabatic perturbations with a nearly scale- 
invariant spectrum. However there are other possibil- 
ities, for example if multiple scalar fields were present 
during inflation, or if the slow-roUedness of the inflaton 
field was temporarily interrupted due to a feature in the 
inflaton potential, etc. Such models can predict power 
spectra that are different from scale- invariance/power- 

law iillMllllllilllSliaiiai). 

Since the details of inflation are unclear, we could use 
precision observations of the CMB and LSS to constrain 
not only the cosmological parameters but also the shape 
of the primordial power spectrum (Pin(fc)) in a non- 
parametric way. The motivation is to test the assump- 
tions that are usually made regarding the shape of Pin{k), 
and to check if there are any features in Pm{k) that can 
thus be unravelled. The aim is to try to see if, as the 
data continue to improve, something new can be learnt 
about interesting cosmologically relevant functions such 
as Pin (A) besides deriving stronger constraints on cosmo- 
logical parameters based on the same set of assumptions. 
We presented an efficient method for reconstructing the 
shape of Pi„(fc) in Mukherjee & Wang (2003c). In this 
paper, we demonstrate that this method is able to ac- 
curately reconstruct an input Pin(fc) by applying it to 
simulated data. 



That Pin(fc) be probed model- independently was ini- 
tially suggested by fsT) and (32), who used conventional 
top-hat binning and linear interpolation respectively in 
the parametrization of Pin(fc). Follow up analysis using 
pre-WMAP data and WMAP data was done by Mukher- 
jee & Wang (2003a) and (H^). In Mukherjee & Wang 
(2003c) we proposed the wavelet expansion method and 
used it on WMAP data. This method is superior to bin- 
ning methods in that it has less correlated errors. The 
method also has various other advantages: due to the si- 
multaneous and adaptive localization property of orthog- 
onal wavelets in both scale and position spaces, the in- 
formation content of most signals is sparsely represented 
in wavelet space with useful information being encoded 
only in a small number of coefficients that are large, while 
the others are small, or in the case of noisy data, consis- 
tent with zero. One can then hope to recover the large 
coefficients more easily and thus constrain the function 
more efficiently. Function reconstruction and denoising 
are thus efficiently achieved with wavelets, as we demon- 
strate further in this paper. 

Since complex projection effects and other degenera- 
cies limit our ability the constrain a function in k space 
using data in multipole space (Tegmark & Zaldarriaga 
2002), the use of a different basis set, that is at the same 
time adaptive to the signal at hand, is certainly a useful 
complementary method to any other model-independent 
method. Since the release of the WMAP first year data 
(Spergel et al. 2003, Peiris et al. 2003), a number of 
analyses with emphasis on the primordial power spec- 
trum have been undertaken (Bridle et al. (2003), Barger, 
Lee & Marfatia (2003), Seljak, McDonald & Makarov 
(2003)). Model-independent methods have been pursued 
in Bridle et al. (2003), (HI El H^; "2^; "2?; © • 

In this paper, we use simulated data to explicitly 
demonstrate the potential and reliability of the wavelet 
expansion method in primordial power spectrum recon- 
struction. Sec. 2 describes our method. Sec. 3 contains 
results. We conclude in Sec. 4. 

2. METHOD 
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CMB data are sensitive to Pin(fc) as follows: 

Q{{b,,l}) = {47t)^J fp,^{k)\ATi{k,T = To)\' (1) 

where the cosmological model dependent transfer func- 
tion AT;(fc, T — To) is an integral over conformal time r 
of the sources which generate CMB temperature fluctu- 
ations, To being the conformal time today. 

As described in l|2ll) we parametrize Pin(fc) by the co- 
efficients of its wavelet expansion. 

J-l 2^-1 

Pinih) = J2T. hi^jAh), (2) 

j=0 1=0 

where ipj^i are the wavelet basis functions, constructed 
from the dilations and translations of a mother function 
The resulting wavelet bases are discrete, compactly 
supported and orthogonal with respect to both the scale 
j and the position I indices, and hence their coefhcients 
provide a complete and non-redundant (hence invertible) 
representation of the function Pin{k). The scale index j 
increases from to J — 1, and wavelet coefficients with 
increasing j represent structure in Pin(fc) on increasingly 
smaller scales, with each scale a factor of 2 finer than the 
previous one. The index Z, which runs from to 2-' — 1 
for each j, denotes the position of the wavelet basis 
within the jth scale. (Note that the properties of the 
basis functions are completely different from the basis 
functions in binning methods). 

We take the sample points ki to be equally spaced in 
log(A:), and use 16 wavelet coefficients to represent Pin(fc) 
in the range 0.0001 < /fc/(Mpc"^) < 0.1, with the spec- 
trum for k < 0.0001 Mpc~^ and k > 0.1 Mpc^^ set to be 
scale-invariant. We use the Daubachies-4 wavelet ^ here 
(@;|2^. In general different discrete orthogonal wavelets 
will give similar results. 

We obtain constraints on the wavelet coefficients of 
Pin(fc), as well as the usual cosmological parameters, 
from simulated data, using a Markov Chain Monte Carlo 
(MCMC) technique to trace the full posterior probabil- 
ity distribution of all the parameters. We use CAMB 
(appropriately modified) and CosmoMC for our compu- 
tations. We simulate cosmic variance limited CMB tem- 
perature power spectrum data in a flat ACDM universe 
with Ho = 70kms~iMpc-\ flbh"^ = 0.022, n^h^ = 0.12 
and T = 0, out to an Imax of 2000, with some functional 
form for the primordial power spectrum. The forms we 
have considered are given in the next section, each nor- 
malized to 1.8 X 10~^. The samples of coefficients can be 
mapped onto corresponding power spectra and the con- 
fidence limits on the power spectra can be obtained at 
each ki , or equivalently constraints on the wavelet coeffi- 
cients can be inverted to obtain constraints on the power 
spectrum with inter-correlations taken into account. 

3. RESULTS 

We consider three different functional forms for the in- 
put Pin(fc), the scale- invariant, saw-tooth, and Gaussian 
dip models. The scale- invariant is given by 

Pi„(fc) = 1. (3) 

^ The number associated with the wavelet corresponds to the or- 
der of the wavelet; the larger the number the smoother the wavelet 
is and the better it is at representing smooth functions or higher 
order polynomials. 



The Gaussian dip model is given by 



Pin(fc) — I ~ ae 



-[log(k/Q.05)f/P 



(4) 



where a = 0.7 and /? = 0.1. The saw-tooth model is 
given by 



Pinik) = ai, k < ki 

= flAT, k > fc^r 



(5) 



k k k'i — ]^ 
a-i-l + : 



h-i 



ki — \ <C k ^ k'i 



where oi = 1, 02 = 0.65, 03 = 1.35 and — 1, with 
the kiS stretching between 0.005 to 0.1 equally spaced 
in log(k). Note that these spectra have been selected in- 
dependently of the discretization of the recovered power 
spectra. What the method efficiently fits for is the cor- 
respondingly discretized input spectra. 

The orthogonality of the wavelet basis functions in k 
space is degraded due to the non-linear model-dependent 
mapping between k and I spaces and due to degeneracies 
with cosmological parameters. As a result it is impor- 
tant to take into account the correlations between the 
constrained wavelet coefficients in deducing constraints 
on the primordial power spectrum. The correlations can 
be obtained with a sufficient number of MCMC samples. 
We obtained about 10^ samples for each case. 

The Ci spectra corresponding to these primordial 
power spectra are shown in Figure 1. Corresponding 
Ci observations are then simulated for the WMAP and 
Planck satellite missions. For the WMAP experiment 
three frequency channels are used with beam FWHM of 
31.8, 21.0 and 13.8 arcmins, and pixel noise of 19.8, 30.0, 
45.6 /xK, respectively. For the Planck experiment, three 
frequency channels with beam FWHM of 10.7, 8.0, and 
5.5 arcmin, and pixel noise of 4.6, 5.4 and 11.7 /iK, re- 
spectively, are used. For both experiments a sky fraction 
of 0.8 is assumed. Instrumental noise of these specifi- 
cations together with cosmic variance uncertainties are 
used to simulate data (see for example Tegmark et al. 
2000). 

Figure 2 shows the input P\n{k) (dashed curve), and 
mean power spectrum reconstructed from simulated Ci 
data together with 2a uncertainties (solid curves). The 
dot-dashed curves give the 5(7 contours. This is shown for 
the scale-invariant, saw-tooth and Gaussian dip models 
of Pin(fc), for WMAP and Planck simulated observations 
of the CMB temperature anisotropy power spectrum. In 
each case the dashed curves shown are discretized ver- 
sions of the corresponding continuous spectra. Figure 1 
shows that the input Pin(fc) is accurately recovered by 
our method. The features are detected at high signifi- 
cance. 

We have used CAMB with its original discretization in 
k and / spaces with interpolation in between; this grid- 
ding can be made smoother to improve the accuracy of 
reconstructions if needed. 

When applied on future data, for example from the 
4-year WMAP or Planck experiments, complemented 
by LSS power spectrum data, we recommend that our 
method be used as follows: 

(a) Follow the procedure described in this paper to de- 
duce constraints on Pin(fc). 

(b) Look at constraints obtained from the larger (more 
significantly constrained) wavelet coefficients, i.e. those 
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that deviate from zero at greater than say Icr, or 2a. 
If tentative evidence of a feature is found in the thus 
denoised Pin(fc) spectrum, i.e. if any wavelet coefficients 
other than the two that together indicate/reconstruct the 
amplitude of a scale-invariant function are found to be 
significantly non-zero, then there is indication of a fea- 
ture or a deviation from scale-invariance. 
(c) If step (b) finds the indication of a feature, one can 
use a parametric form for the kind of deviation from 
scale-invariance that is indicated, and repeat the likeli- 
hood analysis to better fit the feature and to obtain con- 
straints on the relevant parameters. On the theoretical 
side, hopefully the feature will indicate something that 
is theoretically viable in certain scenarios so that it can 
be thus modelled. The relevant parametrization to use 
can also be derived in this way. If evidence calculations 
indicate that scale-invariance or power-law forms remain 
the best explanation for the data we will have achieved 
the task of checking for this model-independently and 
ruling out or constraining more complex physics within 
inflationary models. 

4. CONCLUSIONS 

It is important to reconstruct the primordial power 
spectrum as a free function from cosmological data, the 



aim being to test the assumptions that are usually made 
about Pi„(fc), and to constrain it model- independently 
so as to be able to detect any features in it that may be 
signatures of new physics during inflation. Represent- 
ing Pin(fc) by coefficients of its wavelet expansion is an 
effective and powerful way of approaching this problem 
as described above. We have demonstrated that input 
Pin(fc) are accurately recovered with this method. We 
have presented a prescription for how this method should 
be used on future data. 

We have used CMB temperature data simulated ac- 
cording to WMAP and Planck specifications, with vari- 
ous forms for the primordial power spectra. Our results 
(see Fig. 2) indicate that future data hold great promise 
for probing the detailed shape of the primordial power 
spectrum. 
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Fig. 1. — The plot shows the temperature anisotropy C; spectrum resulting from a scale-invariant (solid curve), saw-tooth (dashed surve) 
and Gaussian dip (dot dashed curve) primordial power spectra. 




Fig. 2. — Each plot shows the input primordial power spectrum (dashed curve), and mean power spectrum reconstructed from simulated 
C; data together with 2<t uncertainties (solid curves) and the 5(7 confidence contours (dot-dashed curves). These are shown for WMAP 
simulated observations in the (a) scale-invariant, (b) saw-tooth and (c) Gaussian dip cases, and similarly for Planck simulated observations 
in (d), (e) and (f). 



